The role of final-state interaction in tensor polarization effects of the reaction \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma d \rightarrow p n \pi ^{0}$$\end{document}γd→pnπ0

Tensor analyzing-power components \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_{20}$$\end{document}T20, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_{21}$$\end{document}T21, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_{22}$$\end{document}T22 for the reaction \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma d\rightarrow np\pi ^0$$\end{document}γd→npπ0 have been studied for the first time in the photon energy range from 280 to 500 MeV. The data are extracted from the experimental statistics accumulated at the VEPP-3 storage ring in 2002–2003. The measured asymmetries are compared with the results of statistical simulations performed with the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma d\rightarrow np\pi ^0$$\end{document}γd→npπ0 amplitude from a spectator model, taking into account corrections for the final-state interaction. The comparison demonstrates quite good agreement between the experimental results and the theory.

Photoproduction of π mesons on nucleons and nuclei is one of the main sources of information about nucleon resonances. The special role of these processes in meson-nuclear physics is due to certain advantages of using photons as sensitive probes. First, the electromagnetic interaction is well understood within the framework of quantum electrodynamics. Second, photons can penetrate deep into nucleons and nuclei, thus making it possible to obtain more complete information about their internal structure. This property distinguishes photoproduction from, for example, scattering of pions, which experience intensive absorption in a nuclear environment due to strong coupling to inelastic channels.
In the region of the photon energies below 500 MeV, several theoretical models were developed to study pion photoproduction on a deuteron, where the impulse approximation is typically used. In this model, the deuteron is actually considered as a system of two nucleons on which the pion is produced like on free nucleons apart from kinematic and binding corrections. The reason for this is the weak binding of the deuteron and its relatively large size. The reaction amplitude is then expressed in terms of photoproduction on a single free nucleon, whereas the second nucleon acts as a spectator. The final-state interaction effects associated with rescattering of the final particles are described in terms of two-body, NN and πN , t matrices.
In most cases, models developed according to this scheme provide consistent description of unpolarized differential and total cross sections, but demonstrate larger deviation for polarization observables. In this regard, polarization characteristics are often used as a sensitive test of various model approaches. It is also well known that polarization measurements give more complete information about the dynamics of the process under study, compared to what can be achieved with an unpolarized cross section. For these reasons, despite the general technical difficulties in conducting experiments with a polarized beam and/or target, measurement of polarization observables are among the most important parts of many research programs aimed at studying photonuclear reactions.
From the theoretical point of view, the influence of the final-state interaction (FSI) has been discussed in detail in a fairly large number of publications [1][2][3][4][5] . It should be noted that, unlike, for example, deuteron photodisintegration γ d → np , where two-body mechanisms are of crucial importance, the incoherent reactions γ d → NNπ are dominated by the single-nucleon mechanism. Although the corrections due to final-state interaction are, www.nature.com/scientificreports/ apparently, most important to the spectator model, their contribution is usually at the level of a few percent of the total cross sections. In particular, as it has been shown in the works cited above, interaction between the final nucleons in the neutral channel npπ 0 leads to an approximately 15 % decrease of the total cross section in the �(1232) region, which is in general agreement with the experimental results 6 . Other mechanisms in which both nucleons are involved (like, e.g., meson-exchange currents) play a minor role, unless the leading mechanisms are suppressed. Despite many theoretical analyses and quite high precision of the available experimental results for γ d → npπ 0 , so far there are no experimental data that could be used to study those FSI features which are directly related to the dynamical properties of the interaction between the final particles. The reason is that the noticeable FSI effects which can be distinguished by comparing theoretical predictions with experimental data are mainly an artifact generated by the plane-wave approximation. For example, due to FSI, a sizable decrease of the pion angular distribution dσ/d� π in the extreme forward direction 1 is simply a trivial consequence of the fact that the resulting cross section contains non-physical contribution from the coherent channel dπ 0 if the plane-wave approximation is used for the final np system. The latter appears due to nonorthogonality of the plane wave e iqr and the wave function φ(r ) of the coupled np system (the deuteron). As demonstrated in 3 , after eliminating this ghost admixture, the remaining FSI effect turns out to be relatively small. In other words, the significant influence of FSI in the reaction γ d → npπ 0 is basically just an unavoidable consequence of using the plane-wave impulse approximation, so that it does not provide any interesting information about the role of np rescattering in this process.
One of the ways to minimize the influence of such ghost FSI effects is to study polarization asymmetries. The latter are expressed in terms of the ratio of different hermitian combinations of spin amplitudes, so that these undesirable effects, entering the numerator and denominator with approximately equal weights, are to a large extent cancelled. In addition, it is obvious that FSI should play a primary role in the kinematic regions that are characterized by a large momentum transfer and where the mechanisms with the participation of both nucleons become important. However, the available experimental data mainly cover the region of small momentum transfer, where the FSI effects are quite insignificant (after eliminating the mentioned influence of non-orthogonality of the wave functions). The only exceptions are the data for γ d → π + nn 7 , obtained near the threshold, and the data for the distribution d 2 σ/(d� π dE nn ) in the same reaction 8 , demonstrating a clear maximum at E nn → 0 coming from the 1 S 0 virtual nn state.
In this work, we present experimental results for three components T 20 , T 21 , and T 22 of the tensor analyzing power for the reaction γ d → pnπ 0 . The data are extracted from the experimental statistics accumulated on the VEPP-3 electron storage ring in 2002-2003. VEPP-3 is an accelerator-storage complex located at the Budker Institute of Nuclear Physics, Novosibirsk. It is designed to accumulate and accelerate electrons and positrons. Presently, VEPP-3 is mostly used in various nuclear physics experiments with internal gas targets and as injector for the VEPP-4 accelerator. It contains the internal target equipment whose key element is the Atomic Beam Source with superconducting sextuple magnets, providing a flux of polarized deuterium atoms with high degree of tensor polarization and a negligibly small vector polarization admixture.
The results of measurements are compared with the theoretical predictions provided by statistical simulation based on the spectator model, which also takes into account the contribution of NN and πN interaction in the final-state.
The paper is organized as follows. In the next section, the method and the formalism used to obtain the components T 2M are described. In section "Results", the data obtained in the present experiment are compared with the results of statistical simulation. A brief discussion of the results and conclusion are given in the last two sections.

Research method
The present experiment was performed with an internal target filled with gaseous deuterium, only for which a high degree of tensor polarization can be obtained. A relatively small thickness of the target was compensated by a high beam current inside the accelerator chamber. A jet of polarized deuterium atoms entered the internal target from an atomic beam source (ABS) installed in the median plane of VEPP-3. At the exit of ABS, the degree of deuteron polarization was close to 100% . A detailed description of the atomic beam source can be found in Ref. 9 .
A general expression for the cross section of pion photoproduction on a tensor-polarized deuteron reads where dσ 0 is the corresponding unpolarized cross section, T 20 , T 21 , T 22 are the components of tensor analyzing power, and d J M ′ M are the Wigner d matrices: The orientation of the target polarization axis with respect to the direction of the photon beam is specified by the polar and azimuthal angles θ H and φ H . The target polarization is determined by the degree of tensor polarization P zz where n 0 is relative population of the deuteron state having zero spin projection on the target-polarization axis. www.nature.com/scientificreports/ The presented work reports on the analysis of data obtained from the experiment conducted in 2002-2003. The recoil proton and neutron were detected by coincidence of the two, lower and upper, systems of detectors (Fig. 1). The lower system, which was used to detect recoil protons, consisted of a set of drift chambers and three plastic scintillators. Recoil neutrons were detected in the upper system by using the time-of-flight method. Six thick scintillators were installed at a distance of 3 m from the target, and a thin scintillation counter was installed at a distance of 1.5 m. The polar angle for recoil protons and neutrons varied between 50 • and 90 • , with their azimuthal angles being within ±30 • for the lower arm and ±12 • for the upper arm.
During the data collection, the polar angle θ H was periodically changed to be one of the three values 54.7 • , 125.3 • and 180 • , whereas the azimuthal angle φ H remained close to 0 • in all cases. The sign of the tensor polarization was switched every 30 s. Such a procedure allowed simultaneous measurement of three asymmetries with respect to the sign reversal of the tensor polarization: is the number of events detected for i-th value ( i = 1, 2, 3 ) of the angle θ H and the target polarization P + zz (P − zz ) . From Eqs. (1) and (4), the required expressions for all three components of the tensor analyzing power T 2M ( M = 0, 1, 2 ) are obtained as The details of the experimental setup, detectors for registering the reaction products, and the methodology adopted for identifying the reaction channel under study are given in Refs. [10][11][12][13][14][15][16] .

Results
The experimental results obtained for all three components T 2M , M = 0, 1, 2 , are presented in Fig. 2 as functions of the laboratory photon energy E γ and of the invariant mass M πn of the π 0 n system. As seen, the asymmetries are quite small and do not exceed 0.2 in absolute value. Because the acceptance corrections for experimentally observed events N + and N − are canceled in the ratio (Eq. 4), they were neglected in the analysis of the experimental data, as well as for the simulated data. The magnitudes of statistical and systematic uncertainties for each data point can be seen in Fig. 2, illustrating a strong dominance of statistical uncertainties. The magnitude of the systematic uncertainties dominates by the uncertainty in determining the degree of the deutron tensor polarization. www.nature.com/scientificreports/ The corresponding statistical-simulation results are plotted in the same figure. It was performed using the Monte-Carlo algorithm described in Refs. 17,18 , which makes it relatively easy to take into account complex boundaries of the experimental kinematic domain, as well as the inhomogeneity of the spatial distribution of the deuteron target.
Statistical simulation was carried out in full accordance with the experimental conditions, including the same constraints on the energies and emission angles of the final-state nucleons. To match the experimental target conditions, the components of the deuteron density matrix were simulated with the same probability, 1/6, for www.nature.com/scientificreports/ each of the six polarization states. Similar to the experimental events, the components T 2M were extracted by using the same formulas (4) and (5) and identical intervals for averaging of the kinematic variables. Such an approach allows a direct comparison of the experimental and theoretical results.
To calculate the reaction amplitude that was embedded into the Monte-Carlo algorithm, we used the γ d → NNπ model developed in Ref. 3 . The amplitude is built within the approximation in which the process on the deuteron is reduced to the sum of single-nucleon photoproduction amplitudes. As already mentioned above, for a process like γ d → NNπ , where the deuteron breaks up, such an approach is called as the spectator model (diagram (a) in Fig. 3). The elementary amplitude γ N → Nπ (shown as a circle in Fig. 3) was taken from the MAID2007 model 19 . It contains contributions from the nucleon born terms, t-channel vector-meson exchange, and a set of s-channel baryon resonances. The latter includes all resonances with masses up to 2 GeV that are classified with four stars in the Review of Particle Physics 20 .
Between the two mechanisms generating final-state interaction, the most important in the kinematical region of the present experiment is the nucleon-nucleon rescattering [diagram (b) in Fig. 3]. This is due to both a larger intensity of the NN interaction, compared to πN , and a fairly high kinetic energy of the final-state nucleons. An additional NN scattering effectively fills the missing-energy balance between the fast active nucleon and the spectator.
The πN interaction is less significant here, primarily because the pion is unable to transfer any large amount of kinetic energy between the nucleons due to its small mass.
For the deuteron wave function, as well as for the NN scattering state, the separable version of the Paris potential from 21 was adopted, in which the partial waves were taken into account up to 2S+1 L J = 3 G 3 . The calculation of the diagram with pion rescattering [diagram (c) in Fig. 3] was carried out using the separable model for the πN amplitude from 22 with all πN partial waves up to l = 2.
Note that in the present calculation, the final-state interaction effects were included only up to the first-order terms in the two-body NN and πN , t matrices, neglecting the higher order terms in the corresponding multiple scattering series. The latter could be taken into account, for example, within a three-body πNN model, as it is done in 23 . As shown in 23 , despite the smallness of the contribution from higher-order multiple-scattering diagrams to the unpolarized cross section, it could still be important. The main reason for neglecting the higherorder terms is in a significant increase of the time required to run statistical simulations with three-body calculations. Thus the question about importance of the terms beyond the first order in the two-body NN and πN rescatterings remains open.
This work also did not address the question about how much the results of simulation in Fig. 3 depend on a model used to construct the elementary amplitude γ N → Nπ . In general, such dependence should not be very crucial because, in the energy range under the consideration, E γ < 500 MeV, the existing multipole analyses for γ N → Nπ differ little from each other. However, because of quite specific kinematic conditions of the experiment (in particular, a large momentum transfer to the spectator nucleon), theoretical results may even be sensitive to small differences between the model amplitudes.
As shown in Fig. 2, taking in to account the final-state interaction effects significantly improves the agreement between the model predictions and the experimental data points, even if some features of the observed tensor asymmetries are not fully reproduced. This is especially important after taking into account quite high sensitivity of T 2M to various model details, as well as to the special kinematic conditions of the present experiment that were discussed in the text above. Such an observation can be reviewed as an indirect confirmation of the general assumption that the spectator model, including NN and πN rescatterings as the next order terms, may be considered as an adequate theory of the process under study.

Conclusion
The present work reports on the first measurements of the tensor analyzing-power components T 2M , M = 0, 1, 2 , for incoherent π 0 photoproduction on a deuteron in the range of the incident-photon energies from 280 to 500 MeV. The experimental results were obtained from analysis of the data accumulated on the VEPP-3 electron storage ring at the Budker Institute for Nuclear Physics in 2002-2003. The present results are compared to the predictions from statistical simulation performed with and without final-state rescattering, which demonstrates that taking such interactions into account significantly improves their agreement.
The results presented in Fig. 2 seems so far to be the only case when the importance of including interaction effects in the incoherent photoproduction of π mesons on the deuteron is unambiguously and quantitatively